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Abstract. - The run-and-tumble dynamics of bacteria, as exhibited by E. coli, offers a simple experimental 
realization of non-Brownian, yet diffusive, particles. Here we present some analytic and numerical results for 
models of the ideal (low-density) limit in which the particles have no hydrodynamic or other interactions and 
hence undergo independent motions. We address three cases: sedimentation under gravity; confinement by a 
harmonic external potential; and rectification by a strip of 'funnel gates' which we model by a zone in which 
tumble rate depends on swim direction. We compare our results with recent experimental and simulation 
literature and highlight similarities and differences with the diffusive motion of colloidal particles. 
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^ Recent years have seen an upsurge of activity at the interface 
^between physics and biology. In many cases, relatively well- 
+^ established physics tools have been used to address problems 

S of major interest to biologists. In other cases problems that 
( might appear arcane to some (though not all) biologists have 
nonetheless posed interesting new physics questions. For in- 
£h stance, the motion of E. coli and similar bacteria involves a se- 
O ries of roughly straight-line 'runs' punctuated by rapid changes 
o f direction ('tumbles') [1]. This run-and-tumble dynamics 
can, in idealized form, serve as a paradigm for non-Brownian 
r ~' diffusive motion, and be used to explore some central concerns 



of nonequilibrium dynamics, such as the origins of phase sepa- 



^-j- ration in systems without detailed balance [3,4]. 

E. coli is a unicellular organism with external flagellar fila- 
C^ments, which form a bundle at one pole of the cell [1]. When 
these filaments synchronously rotate counterclockwise, the or- 
<^> ganism swims in trajectory that is basically straight (with small 
'^effects of rotational Brownian motion) [5]. A stochastic intra- 
^P, cellular event results in a change of activity [6, 7] so that one or 
>- more filaments start rotating clockwise; the organism then tum- 
"kpbles [8,9]. The tumble is of short duration; once all filaments 
start rotating counterclockwise again, a new run begins. 
C3 Supposing the runs to have fixed direction and speed v, and 
the tumbles to occur instantaneously and randomly at rate a 
(with each tumble fully randomizing the direction of motion), 
it is easily shown that in free space such run-and-tumble or- 
ganisms obey the diffusion equation at large length and time 
scales, with diffusivity D = v 2 /ad where d is the dimension- 
ality. Yet, because bacteria are not Brownian particles close to 
thermal equilibrium, one cannot write the usual Einstein rela- 
tions D = kTfi. Here p is the mobility, which controls the 



mean velocity v = pF of a particle subject to an external 
force F; whenever this force is conservative (F = — Vf7), the 
Einstein relation ensures that, for truly Brownian particles, the 
equilibrium mean density p obeys the Boltzmann distribution, 
p oc cxp[- pU/D] = exp[-C//fcT]. 

As will be seen from two of the examples below (sedimen- 
tation and the harmonic trap), forced run-and-tumble particles 
show more complicated behaviour than do Brownian ones. In- 
deed, even when forces are conservative and a flux-free steady- 
state density is recovered, this is not in general a Boltzmann 
distribution for the potential U, even if D is rescaled at will: 
there is no general 'effective temperature' concept. 

This was discussed briefly in [3] which however focussed on 
the case of interacting particles. In the present work we ex- 
plore the noninteracting limit in more detail, presenting a range 
of results that include an exact steady-state calculation of the 
sedimentation decay length in 3D. (The interactions that we 
omit include hydrodynamics [10], which may bring important 
changes even at modest density.) Oddly, we have found very 
little discussion of external force fields in the previous litera- 
ture on noninteracting run-and-tumble particles [11, 12], much 
of which assumes rotational symmetry of the set of possible 
particle velocities. This assumption is inspired by models of 
bacterial chemotaxis, which for good biological reasons take 
the tumble rate, and not the swim speed, to be a function of 
swim direction. But external forces, such as gravity, break this 
symmetry and thus require separate treatment. 

As well as external force fields, we consider below the effect 
of asymmetric obstacles, such as funnel gates [13], arguing 
that these can be modelled by a zone in which there is effective 
asymmetry in the tumbling rates for particles moving in oppo- 



p-1 



J. Tailleur M. E. Cates 



site directions. We show that this viewpoint provides a semi- 
quantitative account of the rectification experiments of [13]. 

Effective temperature: We first consider noninteracting run- 
and-tumble particles in ID with symmetric tumbling rates a 
and swim speeds v. At times long compared to a -1 and length 
scales large compared to t = v/a (the run-length) the motion 
is diffusive, with diffusivity D = v 2 /a. In the presence of an 
external potential U, the speeds of left- and right-going bacteria 
become = v ± /iVU. The diffusion-drift equation for the 
probability density P becomes [3, 1 1] 



d t P=-V-J 

a \ a J 



If the external field only slightly perturbates the swim speed 
of the bacterium, that is if |/xVE/| <C v, and if it does not in- 
duce large gradients of the velocity field on the scale of the 
run length (| V ^ U -C v\), the current becomes at first order 

J = -^VP - fj,(VU)P and the flux-free steady state then 
obeys 

P oc exph&T]; $ = ^ = £ (2) 

Such ID bacteria thus behave as "hot colloids", with an effec- 
tive temperature 1//3; however this breaks down nonperturba- 
tively [3]. We show below that the same is true for d > 1. 

Sedimentation: An obvious physics question concerns the 
sedimentation equilibrium of run-and-tumble bacteria: what 
is the steady-state probability density P(z) to observe a bac- 
terium of buoyant mass m at a height z above a hard wall, in the 
presence of a downward force F = —mg1 A technical com- 
plication is that the singular forces arising at the wall have to 
be described (or replaced by a rule for modifying nearby trajec- 
tories). This will materially influence the profile in a proximal 
region whose height is set by the run length v/a. For simplic- 
ity we here consider only the distal part of the profile at much 
larger z. 

In ID the problem can be solved by the methods of [3, 11] 
which developed diffusion-drift equations for the density that 
nonetheless exactly recover flux-free [1 1] or general [3] steady 
states. A route more readily generalized to d > 1 instead 
addresses the following equations for steady-state probability 
densities which contain no approximation beyond those of the 
defining model: 

P{z)dz = P(c,z)dz (3) 

c=±l 

a f°° 

P{c,z)dz = - / P(z t (z,c,T))dz t e~ aT dT (4) 
z Jo 

Here P(c, z) is the probability density of finding a bacterium 
at height z that swims either upwards (labelled c — 1) or down- 
wards (c = —1) and P(z) is the total probability density. Eq.(4) 
then states that a particle of type c currently (at t = 0, say) in 
the small height interval (z, z + dz) last underwent a tumble at 
some earlier time. (Note that in ID, half the tumbles leave c 
unchanged, whence the 1/2 factor). If this event occurred be- 
tween times t = —t and t = — (r + dr), the bacterium was 



then in a height range (z i: Zi + dzi) where Zi = z — {vc — vt)t. 
Here vt — —fiF so that vc — vt is the net upward velocity of 
a swimmer of type c. (Note that for constant F only, we also 
have dzi/dz = 1.) Thus the contribution to P(c, z)dz of the 
'age-slice' (r, r+dr) is (a/2)dr xP(zi)dzi xcxp[— ar] where 
the three factors represent respectively the probability of a tum- 
ble event occurring in the relevant time interval; the probability 
of finding a bacterium at this time in the required height range 
to later arrive at z; and the survival probability against further 
tumbles at intervening times. Integrating this over r recovers 
(4). 

Summing (4) over c yields a linear integral equation for P(z) 
which can be solved by Fourier transform. Defining P(u>) = 
X^oo P( z ) Gxp(—iujz)dz, (3) and (4) imply 



P(oj) u) [uj(v 2 — v T ) — iav T ] = 



(5) 



P( W )=4%)+Bi(,-^) (6) 

0, one has 



Thanks to the flux-free boundary condition at z 
A = and inverting the Fourier transform yields 

P(z) = p oe - KZ 

VtOL 



(7) 
(8) 



where we have assumed v > vt- In the perturbative regime 
v T <C v, one recovers the result (2). More generally, however, 
although Eq. (7) has the Perrin form (a Boltzmann distribu- 
tion under gravity) Eq. (8) contradicts the result k = vt/D for 
Brownian particles. In particular, the sedimentation length kT x 
tends to zero as vt is raised towards v, while for vt > v both 
species of particle have a net downward motion and nothing 
(short of true Brownian motion, if present) can save the system 
from complete gravitational collapse [3]. 

Unlike the differential equation methods of [3, 1 1], the above 
integral-equation approach generalizes, without further approx- 
imation, to the 3D case. By identical reasoning one finds an 
exact equation for sedimentation of ideal swimmers 



P(z)dz = 2tt J dcP(c,z)dz 



(9) 



P(c, z)dz 



47T 



P(z i (z,c,T)dz I e- aT dT (10) 



where c = cos 9 is the angle between the propulsion direction 
and the vertical, so that the net upward velocity is vc—Vt as be- 
fore, and the definition of Zi(z, c, r) is likewise unchanged. In- 
tegrating (4) over c and proceeding as before via Fourier trans- 
form, one gets for P(uj) = P(z) exp(iojz) 



p{w) = PH 



a + iuiv T \ sin(u) 



/ a + iuv T \ 
exp I du (11) 

U>V Jq V LOV I u 



The r.h.s can be computed, for instance by power expansion of 
the sin function, to yield 



LJV 



arctan 



\a + VjJVt / 







(12) 
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Figure 1: Simulations of 3D run-and-tumble bacteria under gravity. 
Black: full numerics; red: exact solution; blue: quasi-Brownian ap- 
proximation, with vertical offsets introduced for clarity. For v T /v = 
0.7, the deviation from the diffusive approximation is apparent. Inset 
For v T /v — 0.1, the diffusive approximation works very well. 

The constant solution P(uo) = S(uj) is once again forbidden by 
the flux-free boundary condition and one is left with P(lo) = 
S(u}—ujo) where uiqv/o. — a,rcta,n(u>ov / (a+iu)ov T )) ■ Inverting 
the Fourier transform again yields the Perrin form (7), where 
k = iui now has to satisfy 

In ( <;T + v) + a \ = 2nv 
\k(vt ~ v) + a J a 

Just as in ID, linearization of (13) yields the quasi-Brownian 
result (k — vt/D with D = v 2 /3a) at leading order in vt, 
whereas beyond this, although the Perrin form (7) is main- 
tained, its decay rate k again diverges when vt — * v and is infi- 
nite for vt > v. In this regime (barring true Brownian motion) 
gravitational collapse occurs just as found in ID, and for the 
same reason. Interestingly, since (13) and (8) differ, run-and- 
tumble diffusion cannot generally be factorized over orthogo- 
nal spatial directions even if, as here, the forcing is uniaxial and 
translationally invariant. 

For a typical bacterium in terrestrial gravity one expects 
vt — 1/^ms -1 whereas v ~ 20^ms _1 [1]. This means that 
the sedimentation problem (at least in the distal region) remains 
close to the quasi-Brownian limit; however the opposite limit 
could easily be reached using a centrifuge. We hope that our 
exact 3D solution will encourage experiments on bacterial sed- 
imentation. We have also tested our prediction by an explicit 
stochastic simulation of an ensemble of run-and-tumble parti- 
cles, recovering satisfactorily the distal solution (Fig. 1). 

Trapping: Another physics question of possible experimen- 
tal relevance concerns trapped bacteria. Consider, e.g., a har- 
monic potential (U (z) — Xz 2 /2/z) for which vt(z) = —Xz. In 
ID, (3,4) still hold, but solving i = vc — Xz gives in (4) 

Xzi(z, c, r) = vc — (vc — Xz) exp[A"r] (14) 

so that dzi/dz = exp[Ar]. The solution found from the 
diffusion-drift formulation [3] is, for \z\ < v/k, 

P(z) = P(0)[l-(Xz/v) 2 ] a/2X - 1 (15) 



with P(z) = for \z\ > v/X. On substitution in (3,4) with use 
of (14), we confirm (15) to be an exact ID result. 

As noted already in [3] it is not of Boltzmann form (P ~ 
exp[— aXz 2 /2v 2 ]), nor can it be made so by any choice of 'ef- 
fective temperature' (any rescaling of D). As discussed earlier, 
it does however approach the quasi-Boltzmann result if simul- 
taneously (kX/v) 2 <C 1 (defining a region where vt <C v) and 
a/A ^ 1 (i.e. the velocity barely changes over a run length). 
This additional requirement ensures that tumbling is frequent 
on the length scale over which vt changes appreciably and that 
the diffusion-drift approximation (1) holds. (Note that our ex- 
act result (15) is not limited to this case.) The vanishing of P 
for \z\ > v/X arises, as in sedimentation, because in this re- 
gion all particles are moving towards the centre of the trap and 
a finite density cannot be maintained in steady state. 

The breakdown of the 'effective temperature' concept is 
much more severe here than for sedimentation (where the Per- 
rin form was maintained, albeit with an altered k). Whenever 
a/2 X is not large the distribution is strongly non-Gaussian; 
for a = 2X it is flat and for a < 2X it is bimodal, with an 
accumulation of bacteria at the outer edges the trap (Fig. 2). 
This phenomenology has a simple intuitive explanation. First, 
let us note that the motion of right and left going bacteria 
amounts to gradient descent in two effective potentials Ul.r = 
±vz/fi + Az 2 /(2/i) whose minima are Zr l = iu/A. The 
average duration of a run is a -1 whereas after reversal of a 
particle at z = zr^l it takes time In 2/ A ~ A -1 to reach the 
centre. For a ^> A, a bacterium tumbles many times before 
it can cross the trap and it thus diffuses within the average, 
quadratic potential (Ur + U£)/2, giving a Gaussian steady 
state, P ~ cxp[— akz 2 /2v 2 ]. On the other hand, when a -C A, 
bacteria descend their potential to Zr or Zl, spend a long time 
there and then tumble. This generates a bimodal distribution 
whose maxima are at zr and Zl- 

Stochastic simulations of run-and-tumble particles in 3D in- 
dicate a qualitatively similar behaviour for a radial harmonic 
trap vt = — Ar (Fig. 3). Although such a trap (with trapping 
radius large compared to a bacterium's diameter) is hard to re- 
alize experimentally [14], the same physical mechanism (with 
maximal particle density at the outer edge) would apply for 
bacteria trapped, say, in an emulsion droplet. It would be in- 
teresting to look experimentally for qualitative deviations from 
quasi-Boltzmann statistics in such cases. However, as we ad- 
dress elsewhere [15] coupling of the particles to a momentum- 
conserving solvent could in some case drastically alter, or even 
destroy, the flux-free steady states predicted here, as could, on 
larger time scales, the birth and death of bacteria [16]. 

Rectification: An intriguing recent experiment has shown 
how to trap run-and-tumble bacteria by means quite different 
from the external force fields addressed above. This involves 
'rectification' of the random swimming motion so that parti- 
cles move preferentially into the trap [13]. Two rectangular 
microfluidic enclosures contain swimming bacteria (effectively 
in 2D) separated by a wall of funnels (See Fig. 4). Starting 
from a uniform density, the system evolves until at long times 
the two compartments have uniform populations of bacteria, 
but with different densities. The same experiment done with 
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Figure 2: ID run-and-tumble bacteria in a quadratic trap (v = 1, a = 
1). The result of the simulation of run-and-tumble dynamics is dis- 
played in the perturbative case (A = 0.01), the flat case (A = 0.5) 
and the strong trapping regime (A = 5). The later shows bimodality 
as predicted. Also plotted are the gaussian approximation in blue (for 
weak trapping) and the exact result in red (for all cases). They over- 
lap perfectly with the numerics. Inset Differences between numerics 
and : in the weak trapping case, the exact result (red) and the gaussian 
approximation (blue); in the flat case, the exact result (magenta); in 
the strong trapping case, the exact result (cyan). The agreement is in 
all cases very good, the errors being roughly two order of magnitudes 
smaller that the densities. 



non-swimming bacteria results in no rectification. 

Our model neglects hydrodynamics interactions, both be- 
tween bacteria and with the walls [10], and simply assumes that 
a bacterium swims along a wall after colliding with it, until its 
next tumble. (Such motion has been observed in this exper- 
iment [13] and may indeed be of hydrodynamic origin [17].) 
The above fully defines a stochastic model suitable for simu- 
lation (Fig. 4). To create an analytical approximation, we rep- 
resent the wall interactions as an additional, non-stochastic, re- 
orientation process ('forced tumbling') whose rate depends on 
the initial swim direction: while the two particles in Fig. 4 start 
symmetrically either side of the plane of the wall, one will cross 
this plane whereas the other will not. A wall of funnels can thus 
be modelled by a narrow strip in which (say) left-to-right tum- 
bles are favoured over right-to-left ones. 

We next show that this simplified model accounts semi- 
quantitatively for the rectification observed in [13]. First con- 
sider a ID system with densities R(x), L(x) of right- and left- 
moving particles of equal speed v but that tumble with different 
rates an^{x). Neglecting noise terms [3] these obey [11, 12] 



d t R 



-v\7R 



OiR 



R 



OLL 



L 



d t L= vVL 



(16) 



R 



The 1/2 factor in front of the tumbling rates formally account 
for the fact that after a tumble, bacteria choose their new direc- 
tion at random; but any bias in outgoing trajectories can, in ID, 
be absorbed into a redefinition of an,L- 

We now divide the system into three regions: for x € 
[0, L — e] U [L + e, 2L], bacteria tumble symmetrically with a 
constant rate an,L = ct a , whereas for x € [L — e, L + e] bacte- 
ria swimming to the left are converted into right-going bacteria 
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Figure 3: 3D run-and-tumble bacteria in a quadratic trap (v — 1, a = 
1). Left axis: weak trapping (A = 0.01). Numerical results in 
black coincide with the Boltzmann-like approximation in blue P(r) ~ 
r 2 exp(— 3Ar 2 /2). Right axis: strong trapping case (A = 2 in red), 
P(r) /r 2 increases near the edge of the trap. The bacteria accumulate 
on the surface of the sphere of radius v/X. 



with an additional rate a c , i.e. aj, = a a + 2a c , cxr = a Q . 
The wall of funnels is thus replaced by a strip of width 2e in 
which the sole effect is to create a bias in the tumble rates. In- 
troducing the total probability density P = R + L and current 
J = v(R — L), Eqs.(16) are equivalent to the continuity equa- 
tion (1) coupled to one for the current evolution [3, 1 1] 



1 



dtJ 



-DVP + VP-J 



(17) 



with a = (a^ + Q!L)/2; V = v(aL — an)/2a; and D = v 2 /a. 

In the left and right compartments, where or — qj,, the sta- 
tionary solution of (1) and (17) is obviously provided by con- 
stant profiles, which we denote pi and p 2 , respectively. In the 
intermediate strip region however, the solution is an exponen- 
tial increase p(x) oc exp [ia c /4 Matching these three solu- 
tions, one finds for the rectification ratio 
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P2 
Pi 



exp 



2a c e 



(18) 



Allowing for the inverse dependence of D on dimensionality, a 
pseudo-2D version of (18) would then yield 



A p2D = exp 



4a r e 



(19) 



In fact, a fully 2D analysis is also possible for a specific form 
of the tumbling bias. Specifically, if we take constant v and a 
tumble rate a + ai.u for a swimmer moving along unit vector 
u, then the density p(x) obeys in flux-free steady state [11] 



p(x) = p(0) exp 



f ai(x')dx' 
Jo 



(20) 



Note that single-valuedness of the integral requires ot 1 = 
— for some scalar field (j>; without this, there is no flux- 
free steady state [11]. For our problem, taking a i to be a 
constant vector normal to the plane of funnels within the strip 
and zero outside it, we obtain A 2 d — exp(2aie/w). Note that 
the chosen angular dependence of the tumbling rate does not 
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Figure 4: Left: The sides of the enclosures in [13] are 400/im long, 
the arms of the funnels are 27/xm long and their apex angle is 7r/3, the 
funnels are separated by gaps 3.8/im wide. Right: After a collision 
with a wall, a particle tends to swim along it. The wall thus induces a 
tumble with fixed outgoing direction. 

account precisely for what happens close to the funnel and we 
should thus not expect a perfect quantitative agreement with the 
experiment. 

We now compare these theoretical predicitions with our di- 
rect stochastic simulations of 2D mn-and-tumble particles in 
the geometry described previously. We simulate a square cham- 
ber 400/im wide similar to the one used in [13] (Fig. 4, with 
each funnel arm 27/im long and the gaps 3.8/im wide. The api- 
cal angle is 7r/3. In the simulations, when a swimmer hit a wall, 
it aligns with it and proceeds in that direction until the next tum- 
ble. This indeed leads to rectification (Fig. 5), but to compare 
quantitatively with our theoretical model, we must measure the 
effective tumbling asymmetry a c orai. 

Within the ID (or pseudo-2D) approach a lower bound for 
a c is provided by the frequency at which particles swimming 
to the left encounter the apices of the funnels. Numerically, 
this can be obtained as a c — ricoi/t^, where one follows a 
very long trajectory and measures both n co i, the total number 
of such apical collisions, and t e L , the time spent swimming left- 
wards within an enclosing strip of size 2e. We then smoothly 
increased e until ea c converges. In the fully 2D analysis, inte- 
grating the tumble rate a Q + at i.u over each hemisphere gives 
( a L,B.) — a ± 2ai/ir with only half these tumbles chang- 
ing the horizontal direction of swimmers. The corresponding 
lower bound is then provided by ct\ = Trn co i / (2t e L ) . The 
comparisions between numerics and theory (with these limit- 
ing parameter values) are presented in Fig. 5. Increasing the 
swim-speed increases the flux of bacteria colliding with the 
funnels, and thus raises both n co i / t e L and the rectification ra- 
tio. The agreement with theory is semiquantitative, and quite 
satisfactory given the simplicity of the models. A more detailed 
analysis of the angular dependence of the tumbling rate could 
probably improve the quantitative agreement, but might not add 
much further physical insight. 

The 'forced-tumbling' picture of rectification admits several 
parameter-free predictions. Specifically, the rectification ratio 
depends only on the local properties of the funnel wall (gap 
size, apical angle, etc.) and not its overall extent or orientation 
relative to the walls of the chambers; the size and shape of these 
chambers are themselves irrelevant so long as they are large 
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Figure 5: Comparision between theory and simulations (setting a — 
Is -1 for the chamber shown in Fig. 4). Top The rectification ratio A 
measured from numerics (blue, top line), for velocities ranging from 
4/im s" 1 to 60/im s~\ is compared with A P 2d (magenta, middle) and 
Aid (cyan, bottom). The experimental value for the ratio is around 
3 [13], and is reproduced here for v ~ 15/im s _1 . In the experi- 
ment, the velocity was measured to be v = 20 ± 5fim s _1 and the 
tumbling rate a D = 0.75 ± 0.25 s _1 [18]. Bottom As predicted by 
our theory, doubling the system width (from L — 400/im in red to 
L = 800/im in blue) does not change the ratio. When doubling the 
height (to H — 800/im in green) and the number of funnels (from 
13 to 26), the curves overlap until the velocity is large enough to re- 
veal finite size effects. An asymetric chamber (L = 600/im with the 
right cavity twice as wide as the left one) also yields overlapping re- 
sults (cyan). Last, bacteria colliding elastically with the walls are not 
rectified (magenta). 



compared to the run length. (This all follows from the path 
invariance of the integral in (20).) Some simple checks of this 
(doubling the chamber width or height or using an asymetric 
chamber) are made in Fig. 5. 

Detailed balance and symmetry breaking: Our stochastic 
simulations differ significantly from those of Wan et al. [19] 
who simulated a similar geometry. First, these authors assumed 
the tumbling of bacteria to occur at regular rather than Poisso- 
nian time intervals. Our drift-diffusion picture shows this to 
be unimportant; this is fortunate since tumbling events for real 
bacteria are thought to be near-Poissonian [1,2]. Second, Wan 
et al. represented each funnel by a force field deriving from a 
crescent-shaped potential barrier. This force field can overcome 
the swimming force, but cannot change its direction, whereas 
the latter is the key element captured by our own 'forced tum- 
bling' description. Omission of this effect effect may explain 
why, to achieve the experimental value of A ~ 3 [13], these 
authors invoked a run-length v /a a not merely larger than the 
funnel size but comparable to the size L of their entire cham- 
ber. (In this limit, almost every run event encounters the wall 
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of funnels whereas in the experiment only runs starting within 
a strip of width v/a a <C L do so.) Our own simulations instead 
achieve the observed A ~ 3 with experimentally plausible pa- 
rameters (v <~ 15/x ms -1 , a a — ls^ 1 ). 

More fundamentally, it has long been known that to produce 
directed motion out of fluctuations (ratchet effect), one needs to 
break both space-inversion and time-reversal symmetries [20]. 
The spatial symmetry breaking is of course provided by the 
asymmetric shape of the funnels themselves; for a wall per- 
forated by symmetric channels, no rectification can arise [13]. 
Time reversal symmetry is more subtle, however. Whereas the 
motor force of bacteria does break this symmetry, at the coarse- 
grained (diffusion-drift) level, detailed balance is perfectly re- 
stored (albeit with respect to a non-Boltzmann distribution) in 
both the sedimentation and trap examples dealt with above. We 
therefore contend that it is actually the collisions of bacteria 
with the funnel walls that break the time-reversal symmetry in 
the sense required to create a 'ratchet' from the funnel geom- 
etry (Fig. 4). To check this, we simulated the same 2D funnel 
system, with full run-and-tumble dynamics, in a case where the 
entrainment of swimmers upon encountering a wall is replaced 
by a simple elastic collision law. Notably, no rectification what- 
ever is observed in this case (Fig. 5). 

Insofar as the rectification problem is subject to Eq.20, (see 
below for exceptions) its steady states are again subject to de- 
tailed balance once the diffusion-drift coarse graining is ap- 
plied. The dynamics of approach to states states can be ad- 
dressed within the diffusive approximation by neglecting J in 
(17). This yields a Fokker-Planck equation whose late-time re- 
laxation is governed by the first excited state, which for large 
systems is a diffusive mode with relaxation time of order of 
r = L 2 /D. For a velocity of 15/zm s _1 in the geometry of 
Galajda et al., this gives a time scale of the order of 10 min- 
utes, which is quite consistent with the exponential conver- 
gence measured in [13]. 

Steady state fluxes: The restriction «i(x) = — encoun- 
tered in connection with (20) is necessary for steady states in 
which the flux vanishes everywhere. It is easy to envisage ex- 
periments that violate this criterion, for instance if walls of uni- 
form funnel density but both inward and outward orientations 
are used to create a closed shape. If the outward-rectifying sec- 
tions comprise only a small part of the perimeter (the exact cri- 
terion may depend on shape) the mean density within the trap 
will still be higher than outside. Yet, if any such sections are 
present, these will carry a steady outward flux, balanced by an 
influx elsewhere. Thus the steady state contains circulating cur- 
rents. This does not seem to have been looked for in [13]. 

Conclusion: We have calculated and discussed in this let- 
ter several steady-state density distributions within an idealized 
model [11] of noninteracting run-and-tumble bacteria; our re- 
sults included a new exact solution for 3D sedimentation. For 
particles in harmonic traps the mechanism causing extremely 
non-gaussian density profiles at large ratios of run-length to 
trap size was elucidated. The phenomenon of rectification of 
run-and-tumble bacteria by funnel walls was addressed using 
an asymmetric tumble rate within both a ID and 2D model, 
giving semiquantitative agreement with direct stochastic simu- 



lations and, using plausible parameters, with the experiments of 
[13]. We found no rectification in the case where each particle 
collides elastically with a funnel wall, rather than swimming 
parallel to the wall until its next tumble. This suggests, at odds 
with [19], that wall interactions (rather than the intrinsic asym- 
metry of bacterial propulsion mechanisms) are the fundamental 
source of time-reversal asymmetry at the coarse-grained level 
relevant to the rectification experiments of [13]. 
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